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£> ' Abstract 

f*"*- ■ Using the Weyl commutation relations over a finite field ¥ q we introduce a family 

of error-correcting quantum stabilizer codes based on a class of symmetric matrices 
| over ¥ q satisfying certain natural conditions. When q = 2 the existence of a rich 

(Nl ■ class of such symmetric matrices is demonstrated by a simple probabilistic argument 

depending on the Chernoff bound for i.i.d symmetric Bernoulli trials. If, in addition, 
these symmetric matrices are assumed to be circulant it is possible to obtain concrete 
examples by a computer program. The quantum codes thus obtained admit elegant 
encoding circuits. 

=5 

1 Introduction 

> ■ 

^ ■ Let A be a finite abelian group with operation denoted by + and identity 0. We identify A 
with the alphabet of symbols transmitted on a classical communication channel. Consider 
the n-fold cartesian product A n of copies of A. Elements of A n are called words of length n. 
A commonly used group is {0, 1} with addition modulo 2. Let A denote the character group 
of A, the multiplicative group of all homomorphisms from A into the multiplicative group 
of complex numbers of modulus unity. For a = (cq, a 2 , . . . , a n ) T 6 A n we define its weight 
u>(a) to be | cij ^ 0}. We say that a subgroup C n of A n is a t-error correcting group code 
if for every non-zero element x = (xi,X2, • • • , x n ) T in C n , w(x) > 2t + 1. In other words, if 
messages transmitted through a noisy channel are encoded into words from C n and during 
transmission of a word errors at the output occur in at most t positions, then the message 
can be decoded without any error. There is a vast literature on the construction of t-error 
correcting group codes and the reader may find an introduction to this subject and pointers 
to literature in [§, 0. 

A broad class of quantum error correcting codes known as stabilizer codes was introduced 
by Gottesman M and Calderbank et al M (also see & [TBI, [133]). To the best of our knowledge, 
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apart from one computer-generated example [13], all quantum error-correcting codes are 
stabilizer codes. Our aim is to give a new description of the theory of error- correcting 
quantum stabilizer codes. First we introduce some definitions. We choose and fix an N- 
dimensional complex Hilbert space H and consider the unit vectors of TC as pure states of a 
finite level quantum system. If A is a finite abelian group with N elements and {e x \ x G A} 
is an orthonormal basis of TC indexed by elements of A we express it in the Dirac notation 
as \x) = e x . If x = (xi, X2, ■ ■ ■ , x n ) T € A n is a word of length n, we write 

|x) = \x x x 2 ■ ■ ■ x n ) = e xi <g> e X2 <g> . . . <g> e Xn 

where the right-hand side is a product vector in the n-fold tensor product 7i® n of n copies of 
7i. Thus, with the chosen orthonormal basis, every word x in A n is translated into a basis 
state |x) of H® n . 

A quantum code is a subspace C n in H® n . Note that a pure state in TC® n described by 
a unit vector in 7i® n has density matrix \ip){%p\. A density matrix p in 7i® n is a non- 
negative operator of unit trace. In quantum probability, a projection operator E in 7i® n is 
interpreted as an event concerning the quantum system and a density matrix p as a state 
of the quantum system. The probability of the event E in the state p is given by TrpE. 
Messages to be transmitted through a quantum channel are encoded into pure states in TC®" . 
When a pure state or equivalently, a density matrix is transmitted the channel 

output is hypothesized to be a state of the form 

i 

where the operators {Li} belong to a linear subspace A of the algebra of all operators on 
H.® . The operators {L{\ may depend on p, but in order to ensure that p is a density matrix 
it is assumed that LIL^) = 1. By the spectral theorem p can be expressed as 



P 



^2pj\ii>j)(ii)j 



where ipj is an orthonormal set in 7i® n and {pj} is a probability distribution with pj > 
for each j. In other words, the output state p is not necessarily pure even though the input 
state is pure. The operators Lj are called error operators and the linear space A from which 
they come is called the error space. 

Suppose there is a finite family {Mj} of operators in 7i® n satisfying the condition 
£V M'jMj = I and for any output state p with -0 in the code C n , 

j hi 

Then we say that the quantum code C n together with the family {Mj} of 'decoding operators' 
corrects any error induced by {Li} from A. In this context we have the following fundamental 
theorem of Knill and Laflamme f| which gives necessary and sufficient conditions for the 
existence of such a family of decoding operators. 
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Theorem 1.1 || Let A be a family of operators in 7i® n and let C n C 7i®" be a quantum 
code with an orthonormal basis ip 2 , ■ ■ ■ , ipd- Then there exists a finite family {Mj} of 
operators in H® n satisfying the conditions: 

(i) £\ M)M 3 = I; and 

(ii) 

j 

if and only if the following condition holds: 

(i/j p \L\L 2 \ijj q ) = 5 Pt gc(Li, L 2 ) for all Li,L 2 G A, 1 < p, q < d, where c(L 1 ,L 2 ) is a scalar 
independent of p and q and 5 p ^ q is 1 if p = q and otherwise. 

Remark 1.2 The proof of the above theorem is constructive and therefore yields the decoding 
operators in terms of A and the basis ipi, . . . , ipd of C n . In this case we say that C n is an 
A- error correcting quantum code. 

Now we specialize the choice of A. Consider all unitary operators in 7i®" of the form 
U — U\ ® U 2 <8> • ■ ■ <2) U n where each [/, is a unitary operator on H and all but t of the t/j's are 
equal to /. Such a U when operating on ip — tpi ® . . . <^ip n G TC® n produces U\ip) which is an 
n-fold tensor product that differs from if) in at most t places. Denote by At the linear span 
of all such unitary operators U. A quantum code C n is called a t- error correcting quantum 
code if C n is an ^-correcting quantum code. 

2 Quantum codes and subgroups of the error group 

Let (A, +) be a finite abelian group with iV elements and identity denoted by 0. Denote by A 
the character group of A and 7i the iV-dimensional Hilbert space L 2 (A) of all complex- valued 
functions on A, spanned by {|a;)}zeA (where the vector \x) denotes the indicator function 
l x of the singleton {x}). Define the unitary operators U a and V x on 7i for every a e A and 
x e A by 

U a \x) — \x + a), Vx\ x ) = x( x )\ x ) 

where x G A. Then 

X (a)U a V x = V x U a V a G A, X e A. 

These are the Weyl commutation relations between the unitary operators representing A by 
translations and A by multiplications. The family of operators {U a V x \ a G Ax G A} is 
irreducible. 

If a G A n then any element x G A n can be identified with an element of A n so that 

n 

X(a) = Yl Xt(oi) Xi e A,aiE A 

i=l 
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where x = (Xi, ■■■,Xn) and a = (a 1>; . . , a n ). Put C/ a = U ai <g> . . . <g> C/ a „ and V x = V Xl <g> 
• • • ® V Xn . Then {f/a^x | a G A™, % G A n } is again an irreducible family of unitary operators 
satisfying the Weyl commutation relations 

x(a)U a V x = V x U a V a G A n , \ G A n . 

In the Hilbert space of all linear operators on 7i®" equipped with the scalar product 
(X\Y) = TrX^Y the set {N~ n / 2 U a V x | a G A n , x £ ^4 n } is an orthonormal basis. The 
weight wt(a, x) of a pair (a, x) G A n x A n is defined to be #{2 | 1 < i < n, (a«, xi) (0, 1)}, 
where a = (ai, 02, ■ ■ ■ , a n ) and x = (xi, • • • , Xn)- The irreducibility of {U a V x | a G A n , x G 
A n } implies that {U a V x | a G A n , x G A™, u>£(a, x) < t} spans A t - The Knill-Laflamme 
theorem for ^-correcting quantum codes assumes the following form which can be readily 
derived from Theorem |1 . 1| . 

Theorem 2.1 C n C L 2 (A)® n is a t-error correcting quantum code if and only if C n has an 
orthonormal basis ipi,' l p2, ■ ■ • , satisfying the following conditions: 
For every (a, x) £ A n x A n such that wt(a, x) < 2t 

(i) (i'i\U a V x \'ip j ) = ifi^j, and 

(ii) (ipilUaV-xlipi) is a scalar independent of ipi for i — 1,2, . . . , d. 

2tyi 

Let / be the least positive integer such that la = for all a G A and let u = e~ . We 
define the error group as the following finite group of unitary operators in L 2 (A)® n . 

£ = {uj l U a V x | < i < I - 1, a G A n , X e A n }. 

The group £ has a natural action on the Hilbert space L 2 (A)® n defined by: 

*7 a |x) = |x + a), V x \x) = x(x)|x). 

Subspaces of L 2 (A)® n that are point- wise fixed by some subgroup of the error group £ are 
called stabilizer codes. 

Let 5 be a subgroup of £. Denote by C(S) the subspace of L 2 (A)® n that is point-wise 
stabilized by S. More precisely, 

C(S) = G L 2 (A)®" \U^ = i> V C/ G 5}. 

Lemma 2.2 C(S) ^ if and only if S is an abelian subgroup of £ such that uj 1 I G" S for 
i 7^ 0. Furthermore, when C(S) ^ the dimension of C(S) is j^A n /#«S. 

Proof. Suppose uj 1 I G S for some « ^ 0. For any ip G C(S) we have a;*/^ = ip , which 
implies ^ = 0. Hence C(S) = 0. 

It follows from the Weyl commutation relations that two elements u/C/aK* and ujiUbVp 
in S commute if and only if ct(b) = /3(a). Now, let ip G C(5). We have 
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Applying the commutation relations we can see that the above equation holds for a ip 7^ 
if and only if a(b) = /3(a). Thus, C(S) 7^ if and only if S is abelian and u l I G" S for i ^ 0. 

Now, let S be an abelian subgroup of £ such that u l I G" S for % ^ 0. Define the projection 
operator 

^ t/es 

Since TrUaVp = unless (a,/3) = (0, 1) it follows that Tr(P) = #A n /#S. It is easy to see 
that P is the projection onto C(S). Thus, the dimension of C(S) is Tr(P) = j^A n jj^S. This 
completes the proof. ■ 



Next, we state Theorem [2.1| in a form that will give the criteria for constructing t-error 



correcting quantum stabilizer codes. Let Z(S) denote the centralizer of S in £, i.e., 

Z(S) = {U e£\UU' = U'U V U' e S}. 

Theorem 2.3 Let S be an abelian subgroup of the error group £ such that u l I is not in 
S for i ^ 0. Then C{S) is a t-error correcting quantum code if wt(a,a>) > 2t for each 

^UaVa G Z{S) \S. 

Proof. Suppose wt(a, a) > 2t for each ^UaYa G Z(S) \ S. Now, by the previous lemma 
C(S) is a subspace of L (A) of dimension j^A n jj^S = d. Let tpi, . . . , ipd be an orthonormal 
basis of C(S). Consider a (a, x) <E A n x A n with the property that wt(a, x) < 2t. We check 



the Knill-Laflamme conditions (Theorem 2A). There are two cases: 

(a) If u l U a V x G S for some % > then 

(ipj\uJ*U a V x \4) k ) = (ipj I tp k ) = S jk , 1 < j, k < d. 
Thus, (ipjlUaY^tpk) — ^~ l fijk\ < j,k < d, where 5j k is the Kronecker delta function. 

(b) If u} i U a V x G" S for each i > 0, then since wt(a, x) < 2t, u i U a V x G" Z(S) for each i > 
by the assumption. Let ip G C(S) and ufV-^f^ be some element of S. Then we can 
write (■0|[/ a V^ : |"0) as (uj r UbVp^\U a V x \uj r Ui } V/3'il)), which can be simplified to get the 
following 



U>\U«V X \1>) = /3(a)x(b)(^|^ x |^). (2) 

Since u/{7 a V^ G" Z(iS>) for each i > 0, for some u r UbVp G 5 we must have /3(a) 7^ x(b). 
This choice of u r Ux,Vp G 5 yields {ip\U a V x \ip) = 0. 



At this point it is useful to introduce a standard notation using which it is convenient to 
describe quantum stabilizer codes. Let S be an abelian subgroup of £ with centralizer Z(S). 
The minimum distance d(S) is defined to be the minimum of 

{wt(a, a) I w'{/ a y a G Z(S) \ S}. 
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When A is the additive abelian group of the finite field ¥ q we define an [[n, k, d]] q quantum 
stabilizer code to be a g fc -dimensional subspace C(S) of L 2 (¥ q )® n , where S is an abelian 
subgroup of £ with d(S) > d and cardinality q n ~ k . 

By Theorem |2.3| it follows that an [[n, k, d]] q quantum stabilizer code is a [(d — l)/2j-error 
correcting quantum code. 

Remark 2.4 Let S be an abelian subgroup of £ such that u l I G" S for every i ^ 0. This is 
equivalent to demanding that S is an abelian subgroup of £ such that for any a G A n and 
X G A n the operator oj l U a V x can be in S for at most one i : < i < I — 1, Thus S has the 
form 

S = {p(ai,x)U a V x \( a ,x)eS} 

where S C A n x A n is a subgroup satisfying x(a') = %'(a) for any (a, x)i ( a 'x') £ S an d P 
is a function on S with values in {u l | < % < I — 1}. 



3 Quantum stabilizer codes in the finite field setting 

In order to construct stabilizer quantum codes, we need to study abelian subgroups S of 
£ such that elements in Z{S) \S have large weight. We choose A to be a finite field ¥ q , 
q = p r for some prime p. In particular, the Hilbert space in which we seek stabilizer codes 
is L 2 (W q )® n . Since ¥ q is an abelian group under its addition operation with each nonzero 
element of order p, it follows that every nontrivial character of ¥ q is of order p. Choose a 
nontrivial character uj G F q . Then every other character lu' G ¥ q is of the form u a where 
^a(x) = Cu(ax) for all x G ¥ q . Likewise, every character in ¥ q is of the form oj a where 
ci; a (x) = u)(a-x) for all x G ¥ q , where a - x is the inner product ^ a^, for a = (a%, . . . , a n ) T 
and x = (xx, . . . , x n ) T . 

If we identify ¥ q with ¥ q , we can index the elements of the error group £ as a/[/ a Vb, 
< i < p — 1, and a, b G F" where Vb now stands for the operator V x with x = u h- Thus, 
£ is rewritten as 

£ = {Ju a V h |0<i<p-l,a,be ¥ q }. 

Notice that £ is a finite group of cardinality pq 2n . The Weyl commutation relations take the 
following form 

u(b-a)U a V h = V h U a V a,beF;'. 

If S is a subgroup of £ it is readily seen that S is abelian if and only if for any two 
elements cu*[/ a Vb and to^UcVd in S we have a ■ d = b ■ c. For (a, b) G F™ x F™, define 
wt(a, b) = | (cij, bi) 7^ (0, 0)}, where a = (a±, . . . , a n ) and b = (bi, . . . , b n ). Let S C 
F™ x F™ be a subgroup for which a • d = b • c for all (a, b), (c, d) G S. Define 

S ±3 = {(a, b) G F" x F™ | a • d - b • c = for all (c, d) G S}. 



Lemma 2.2 and the Knill-Laflamme conditions can be restated as follows. 



Lemma 3.1 Let S C F" x F™ be a subgroup for which a • d = b ■ c for all (a, b), (c, d) G S. 
Suppose p : S — > {uj % \ < i < p— 1} is a function such that S = {p(a, h)U a Vb \ (a, b) G S} 
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is an abelian subgroup of £. Then C(S) C L 2 (F g )® n is a quantum stabilizer code of dimension 
q n /#S. Furthermore, if wt(&,h) > 2t for all nonzero elements (a, b) G S ±s \ S then C(S) 
is a t-error correcting quantum stabilizer code. 

Thus the problem is to find subgroups 5* of ¥ q x ¥ q such that a • d = b ■ c for all 
(a, b),(c,d) G S and wt(a, b) is large for nonzero elements (a, b) G S ±s \ S. The other 
problem is to ensure that we can build an abelian subgroup S of £ by picking a suitable 
p : S — > {uj % | < i < p — 1} such that S = {p(a, b)£/ a Vb | (a, b) G S}. To this end, we 
formulate an approach. 

Let V be an m-dimensional vector space over ¥ q for a positive integer m. Let L : V — > F™ 
and M : V — > F™ be two linear transformations. Thus, L and M can be written as n x m 
matrices over F ? . We restrict attention to abelian subgroups of £ that are of the form 

{p(v)U Lv V M v I v g V}. 

Two elements P('Vi)Ulv 1 Vmv 1 , p(v2)L r L V2 VMv 2 commute precisely when 



v^L T M Vl = vfL T Mv 2 V v 1; v 2 G V (3) 

and 

P(V1 + V 2 ) Trr ., \ ,~, /,„T r T 



p(vi)p(v 2 ) 



uiyiL 1 Mv^ =u(y[L 1 Mv 2 ) V v 1; v 2 GV. (4) 



Equation (||) will hold if we choose L and M such that M T L = L T M (i.e. L and M are 
such that M T L is symmetric). 

Writing p(v) = u)(g(v)), for some function q : V — > ¥ q , Equation (§) assumes the form 

?(vi + v 2 ) - g(v x ) - g(v 2 ) = v^L T Mv 1 = vlL T Mv 2 V v 1; v 2 G V. 

For p ^ 2 we can choose g to be the quadratic form |v T M T Lv. For p = 2 the problem of 
recovering a suitable quadratic form as a solution to the above equation is more difficult. 

For the purpose of this article, we look for special solutions: we demand that L T M be 
expressible as D + D T for some matrix D over ¥ q , which implies that L T M is a symmetric 
matrix with diagonal entries as scalar multiples of 2. For example, we can choose D to be 
an upper diagonal matrix. Then g(v) = v r Dv is a solution to Equation ([|). We summarize 
this below. 

Lemma 3.2 Let V be a finite dimensional vector space, and L : V — > F^ and M : V — >• F™ 
6e £u>o linear transformations such that M T L is symmetric and of the form D + D T for a 
linear map D : V —>■ ¥ q . Then 

S = {0(v T Dv)U Lv V Mv | v G V} 
is an abelian subgroup of the error group £ on L 2 (¥ q )® n . 
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An element u l U x V y of £ is in Z{S) if and only if v T M T x = \ T L T y for all v G V. 
Equivalently, u l U^V y G #(«S) if and only if M T x = L T y. 

From the Knill-Laflamme conditions as stated in Theorem |2.3| , C(S) is a t-error correcting 
quantum code with S defined as above if for any (x, y) G F™ x ¥ q , the condition M T x = L T y 
implies that either x = Lv and y = Mv for some v G V or wt(pc, y) > It. 

If ¥ q has characteristic different from 2 there is a partial converse to Lemma |3.2| : Suppose 
C(S) is some stabilizer code in L 2 (F g )®" where S = {p(a, b)[/ a Vb | (a, b) G 5} for some 
additive subgroup 5* such that a ■ d = b ■ c for all (a, b), (c, d) G S. Let #S = q r and 
(aibx), (a 2 ,b 2 ) . . . , (a r b r ) be an independent generating set for S. Then S\ = {a G ¥ q | 
3beFJ : (a, b) G S} and 5* 2 = {b G F™ | 3a G F™ : (a, b) G 5} are linear subspaces of 
F£. Let ei, e 2 , . . . , e r be the standard basis for ¥ r q . Define L:W q ^> F™ and M : FJ -»• F™ by 
letting Lej = a^ and Me^ = bj for z = 1, 2, . . . , r. Since , a - d = b ■ c for all (a, b), (c, d) G S, 
it follows that L T M is symmetric. Suppose F 9 is of characteristic p ^ 2. For (a, b) G S*, 
let v G FJ be such that Lv = a and Mv = b and define p'(v) = |v T L T Mv. It is easy to 
check that p(a, b) = p'(v) + c ■ v, for some c G ¥ q . More precisely, we have the following 
proposition. 

Proposition 3.3 Suppose ¥ q has characteristic different from 2 andC(S) is some stabilizer 
code in L 2 (¥ q )®" of dimension q n ~ r . Then there are linear transformations L : F^ — > ¥ q and 
M : ¥ q — > ¥™ such that L T M is symmetric and there is a c G ¥ q such that 
S = {cD(|v T L T Mv + c ■ v)U Lv V Mv |vGF;}. 

We can derive the following proposition from Lemma |3.2| . 

Proposition 3.4 Let L : F™ -1 — > ¥ q be an injective linear map with range C = {a G F" | 
SiLi a i = 0} an d M = M'L for some symmetric linear map M' : ¥ q — ► F™ of the form 
M' = D + D T . Then 

(i) S = {cD(a T Da.)U a VM'a \ a E C} is an abelian subgroup of £ . 

(ii) C(S) is t-error correcting if for any (x, y) G F" x F™ the condition y — Af'x G Co- 
mplies that either x G C and y = M'x or u>£(x, y) > 2t. 

Let S = {u(sl t Da)U a VLa aGC}, where L = D + D T , L and D are n x n matrices 
over ¥ q , and C is a subspace of F™. As already observed S is an abelian subgroup of £. Our 
next goal is to give an orthonormal basis for C{S). Notice that C(S) is a g^/^C-dimensional 
subspace of L 2 (F 9 )®". Since C is an additive subgroup of F™, it suggests that an orthonormal 
basis for C(S) can be indexed by the cosets of C in ¥ q . It suffices to describe unit vectors 
l^c+x) £ L 2 (W q )® n that have disjoint support in W q , and show that each \ipc+x) is fixed by 
S, where x runs over a set of distinct coset representatives of C in ¥ q . Define 

|^c + x) = -^^^(a T J Da)cu(a T Lx)|a + x) (5) 

for each coset C + x as x runs over a set of distinct coset representatives of C in F™. The 
vectors |V>c+x) have unit norm, and as they have mutually disjoint supports, they form an 
orthonormal set of q n jj^C vectors in L 2 (F g )® n . It can be easily verified that S fixes each 
l^c+x)- We summarize or observations below. 
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Proposition 3.5 Let S = {u(a T Da)U a VL a | a G C}, where L = D + D T , L and D are 

nx n matrices over¥ q , and C is a subspace o/F™. Then the collection of vectors {\ipc+x.}} 
defined as 

iVto+x) = -™ 5 Y]a)(a T Da)cj(a T Lx)|a + x) 

/or eac/i coset C + x as x runs over a set of distinct coset representatives of C in F™ ; is an 
orthonormal basis for C(S) . In particular, dimC(S) = q n - dimC . 

Remark 3.6 C(S) is an [[n, k, d]] q quantum code if it has dimension q k and d(S) > d. In 
line with classical coding theory we can define the rate of an [[n, k, d]] q quantum code as k/n 
and relative distance as d/n. It is clearly desirable to design quantum codes with large rates 
and relative distance. An [[n, k, d]] q quantum code C(S) is a pure code if the corresponding 
centralizer subgroup Z(S) has the property that wt(&,h) > d for each k/?7 a Vb G Z(S), 
(a, b) ^ (0,0). (Notice that this is a stronger property than guaranteed by Theorem \2.$ ). 
Given a pure quantum stabilizer code, the following simple method can be used for deriving 
new quantum codes. 

Suppose we have an [[n, k, d]] q pure quantum code, with a small k and large d. From such 
a code we can construct an [[n — 1, k + 1, d — l]] g quantum code that is again pure, by the 
technique of puncturing S to yield an additive subgroup S' o/F™ -1 x F™ _1 ; of size still q n ~ k 
and distance at least d{S) — 1. The idea of punctured classical codes (see McWilliams and 
Sloane j^j) can be adapted to punctured pure quantum stabilizer codes following j^j where it 
is shown for q = 2. A repeated application of puncturing will give [[n — k',k + k', d — k']] q 
codes for different choices of k' . 



4 A class of stabilizer codes 

First choose and fix the following subspace C of F": 

C = {(a 1 ,..,a n ) T 6FJ|^a 1 =0}. 

i 

The subspace C is invariant under the cyclic shift permutation a : i t— > (i + l)mod n. 
Thus, C 1 - = {(a, . . . , a) G F™ | a G F^} is also invariant under a. An n x n matrix L over F 9 
is said to be circulant if for % = 2, . . . , n, the i th row of L is obtained by applying a l ~ l to the 
first row. 

Let S = {cD(a T Da)f/ a Vi a | a G C}, where L = D + D T is an n x n matrix over ¥ q and 
C is as chosen above. We further specialize our construction by choosing L to be an n x n 
symmetric circulant matrix with entries from {0, 1} and with all diagonal entries 0. Let 
ei = (1, 0, . . . , 0) G F™. For such an L observe that a T Lex = a T D T ei for a G C. Then the 
orthonormal basis {|V>c+x)} for C(S) (as described in Equation (§)) can be written in the 
following form: 

|^c+ cei ) = -L^ Vw((a T + ce[) J D(a + ce 1 ))|a + ce 1 ), c G ¥ q . (6) 
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In particular, for q = 2 the above stabilizer code has a neat encoding circuit that we describe 
in Figure ?? in the appendix. 

As an example of stabilizer codes given by Equation (||), we now describe a [[5, 1,3]] 9 
quantum code for every finite field ¥ q . In particular, for q — 2, the [[5, 1,3]]2 code is the 
Laflamme code which was originally obtained by a computer search |J. Let L5 be the 
following symmetric circulant matrix in F^ x5 . 

/ 1 1 \ 

1 1 

1 1 
110 

\ 1 1 J 

and C = {(01, . . . , as) T G F^ | ^\ a» = 0}. It can be checked that S — {(a, L 5 sl) a G C} 
is an additive subgroup of F^ x F^ such that d(S) > 3. Thus, by Theorem ^73] , C(S) is a 
[[5, 1, 3]] 9 quantum code for every finite field ¥ q . The encoding circuit for the [[5, 1, 3]] 2 can 
be obtained easily from the general encoding circuit already described for codes given by 
Equation (|6]). 

For a vector c G F£, let ac G F£ denote the vector obtained by a cyclic shift of c. An 
n x n circulant matrix with first column c eFj can be conveiently written as 

( c ere ... cr n_1 c ) 

We give two more examples of quantum codes defined using circulant matrices. 

First, there is a [[13, 1, 5]]2 quantum code defined by a 13 x 13 circulant matrix L13 over 
F2, whose first column is 

c = (0,0,1,1,0,0,0,0,0,0,1,1,0) T . 

As defined, C = { (a l5 . . . , a 13 ) T G Fg 3 | J2i a i = 0}- ^ can checked (with the help of a 
computer program) that S = {(a, L 13 a) | a G C} is an additive subgroup of F^ 3 x ¥ g 3 such 
that d(S) > 5. Thus, by Theorem |2.3|, C(S) is a pure [[13, 1,5]]2 quantum code. 

Similarly, there is a [[21, 1, 7}}2 quantum code defined by a 21 x 21 circulant matrix L21 
over F2, whose first column is 

c = (0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 1, If 

As before, C = {(ax, ■ ■ ■ , ci2i) T G Fif | J2i a * = 0}- ^ can be checked using a computer 
program that S = {(a, L 2 ia) | a G C} is an additive subgroup of Fif x F^ 1 such that 
d(S) > 7. Thus, by Theorem |2~3| , C(S) is a pure [[21, 1, 7}} 2 quantum code. 

If k = 1, it is interesting to note that for n = 5, 13 and 21, the best achievable minimum 
distance is = 3, 5, and 7 respectively. 

5 Existence of good stabilizer codes 

Using a probabilistic argument we show that there is a number a > and a natural number 
n a such that for each n > n a there exists a [[n, 1, L ara J]]2 pure quantum stabilizer code. Now, 
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as observed in Remark |3l| given (3 such that < (3 < a, by the method of punctured codes 
we can obtain a family of [[[(1 — /3)n\, \J3n\, [(a — (3)n\]] 2 quantum codes for all n > n a . 
These are good quantum codes with constant rate (3/(1 — 0) and constant relative distance 
(a-(3)/(l-(3). 

We first recall a particular form of the Chernoff bounds for bounding the probability that 
a random variable deviates far from its expectation. 

Theorem 5.1 || Theorem 4.2, page 70] Let X\,X 2 , ■ ■ ■ ,X n be independent Bernoulli ran- 
dom variables such that for each i, Pr[Aj — 1] — p and Pr[Xj = 0] = 1 — p, for < p < 1. 
Let X = J2i Xi and let \i denote the expectation E[X]. Then for < 5 < 1 



Pr[X < (1 - 6)n] < e 



-fj,6 2 /2 



Our existence proof for stabilizer codes will be guided by Lemma [372 
As before, we first choose and fix the following subspace C of F 

C = {(a 1 ,...,a n )eW^\^2a i = 0}. 



Definition 5.2 An n x n matrix R over F 2 is said to be a-good if the following conditions 
are true. 

(i) The sum of every [cm J columns of R has weight at least an. 
(ii) The sum of every \_an\ rows of R has weight at least an. 
(Hi) The sum of every [an\ columns of R has weight at most (1 — a)n. 
(iv) The sum of every [an\ rows of R has weight at most (1 — a)n. 

As in classical coding theory ||, given a vector a = (a\, a 2 , . . . , a n ) T G F" we denote 
| dj 7^ 0} by w(a). The next proposition describes a way of constructing stabilizer codes 
from good matrices. 

Theorem 5.3 For < a < 1, suppose R is an n x n a-good matrix over F 2 . Let L be the 
following 2n x 2n symmetric matrix over ¥ 2 : 

R 
R T 

// we write L = D + D T , where D is the upper triangular matrix with zeros on the principal 
diagonal, and define the abelian subgroup S of £ as S = {uj(a T D a)U a VL a \ a G C}, then 
C(S) is a [[2n, 1, L Q;n J]]2 pure stabilizer code. 
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Proof. From Lemma 3^2, we know that C(S) is a [[2n, 1, [ an J]]2 stabilizer code if for any 
(x, y) G F| n x Fg™, the condition y — Lx G C 1 - implies that either x G C and y = Lx 
or wt(pc, y) > an. It is easy to check that the assumptions about R in Definition [5.2| , in 
fact, guarantees a stronger property: for any nonzero vector x G F 2 such that u>(x) < an, 
the assumptions (i) and (iii) imply that an < w(Rx) < (1 — a)n. Similarly, assumptions 
(ii) and (iv) imply that an < w(R T x) < (1 — a)n. Putting these together, it follows that 
an < w(Lx) < (1 — a)n if w(x) < an for x G F^™. 

Since C -1 = {(1, 1, ... , 1) T , (0, 0, ... , 0) T }, we can see that the above observation implies 
that wt(x,y) > an if (0,0) ^ (x,y) G ¥ 2 2 n x F 2 n , and y — Lx. G C 1 . It follows that C{S) is 
a [[2n, 1, |_CK7-iJ ]] 2 pure stabilizer code. This completes the proof. ■ 

We now show the existence of n x n matrices over F2 that fulfill the conditions of Theo- 



rem 5.3. 



Lemma 5.4 Let Rij, 1 < i,j < n be independent identically distributed random variables 
taking values in {0, 1} such that Pr[Rij = 1] = 1/2, 1 < i,j < n. Let R be the uniformly 
distributed nx n random matrix over¥ 2 whose ij th entry is the random variable R^. There 
exist constants a > and n a > such that 

Pt[R is a-good ] > 0. 

Proof. Let BAD denote the event that R is not a-good. Let ri, r-2, . . . , r n be the rows of 
R and c 1; c 2 , . . . , c n be the columns of R. For any subset 5" C {1, 2, . . . , n} with 1 < #5 < 
an, we define Es, D s , A s , B s as the events u;(@ ig5 rj) < an, w(@ igS ,r-j) > (1 — a)n, 
w ((&ies c *) < cm ' u '(©ie5 c i) > — a ) n respectively. Then BAD can be written as follows 

BAD = [J A s UB s UDsUE s . (7) 

Sc[n],l<#S<an 

We analyze As for a fixed S. Let © ie5 Cj = (xi, X2, . . . , x n ) T . Since Rij 1 < i, j < n are all 
independent random variables taking values in F 2 , x±, X2, ■ ■ ■ , x n are n independent uniformly 
distributed random variables taking values in F 2 . We will use Chernoff bounds as given in 
Theorem |5.1| to analyze the random variable #{i | Xj = 1,1 < % < n}. Let X = Yli=i x i- 



Then E[X] = n/2. Applying Theorem |5.1| we get 



n( 1 — za ) 

Pr[A s ] = Pr[X < an] < e 2 — . 

Notice that under F 2 addition 1 +xi, 1 + X2, ■ ■ ■ , 1 + x n are also n independent uniformly 
distributed random variables taking values in F 2 . Thus, by Chernoff bounds we again obtain 

n(l— 2q) 2 tc(1-2q) 2 

P r i-BsJ < e 4 Likewise, Pr[i?5j and Pr[D s J are also bounded above by e * 
Putting these together with the definition of BAD in Equation (|7|) we get 

[cmj , . 

n(l-2a) 2 ^ — -v / n\ n(l-2a) 2 „, \ 

Pr[BAD] < 4e"^^ ■ ^ ( . J < 4 e -^^2 nH(Q) 

i=i ^ J ' 
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where H(a) = — a (log a) — (1 — a) log(l — a). To ensure that Pr[BAD] < 1, it suffices to 
pick a < 1/4 such that H(a) < (loge)3/8 — 2/n, which can be done by choosing n larger 
than some constant n a and a > sufficiently small. ■ 



From Theorem 5.3, Lemma 5.4, and Remark p. 6| we can immediately deduce the following. 



Corollary 5.5 There are constants a > and n a > such that for each n > n a there is 
a [[n, 1, L an J]]2 pure quantum stabilizer code. Furthermore, for any [3 such that < (3 < a, 
and n > n a there is a [[[(1 — /3)n\, |_/3nJ, \_(a — /3)nJ]] 2 pure quantum stabilizer code. 



Remark 5.6 The above existence argument can be easily extended to stabilizer codes over 
any finite field. More precisely, for ¥ q there are constants a > and n a > such that for 
each n > n a there is a [[n, 1, [ an J]]g pure quantum stabilizer code. Also, given a (3 such that 
< (3 < a, and n> n a there is a [[[(1 — (3)n\, \_(3n\, [(a — (3)n\]] q pure quantum code. 

Acknowledgment We thank Piyush P Kurur for his help with the computer search for 
the codes given in Section |j. 
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The encoding circuit 



We describe the encoding circuit for the stabilizer code with orthonormal basis defined by 
Equation 6 in the figures below. The elements that we use to build our encoding circuit are 
the Hadamard gate, the C-NOT gate and the Z gate. We first describe these gates. Then 
we give the four main components B, H® n , D, and E and finally, the complete encoding 
circuit in Figure 8. 





H 









(|0) + |l))A/2 



Figure 1: The Hadamard gate. 
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Figure 2: The controlled NOT gate. 



Figure 3: The controlled Z gate. 
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Figure 4: B: computing a + be\. 
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Figure 6: The linear map D. 



Figure 7: Component E. 
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Figure 8: The encoding circuit. 
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